From quantum to classical description of intense laser-atom 
physics with Bohmian trajectories 

X. Y. Lai\ Qingyu Cai^O and M. S. Zhan^'^ 
^ State Key Laboratory of Magnetic Resonances and Atomic and Molecular Physics, 
Wuhan Institute of Physics and Mathematics, 
The Chinese Academy of Sciences, Wuhan 430071, China and 
"^Center for Cold Atom Physics, The Chinese Academy of Sciences, Wuhan 430071, China 

Abstract 

In this paper, Bohmian mechanics is introduced to the intense laser-atom physics. The motion 
of atomic electron in intense laser field is obtained from the Bohm-Newton equation. We find the 
quantum potential that dominates the quantum effect of a physical system becomes negligible as 
the electron is driven far away from the parent ion by the intense laser field, i.e. the behavior of 
the electron smoothly trends to be classical soon after the electron was ionized. Our numerical 
calculations present a direct positive evidence for the semiclassical trajectory methods in the intense 
laser-atom physics where the motion of the ionized electron is treated by the classical mechanics, 
while quantum mechanics is needed before the ionization. 
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I. INTRODUCTION 



In recent years, intense laser-atom physics has received more and more attention 
due to the nonhnear multiphoton phenomena, such as high-order harmonic generation and 
above-threshold ionization. Without doubt, such multiphoton phenomena could be repro- 
duced by resolving the time-dependent Schrodinger equation (TDSE). In order to understand 



such phenomena intuitively, some semic 



arguments, have been proposed 



Qflfl 



assical approaches, mixing classical and quantum 



5j. In the semiclassical approaches, the motion of 
ionized electrons is treated by classical dynamics directly, while quantum mechanics is used 
before the ionization. These semiclassical approaches have obtained much success in intense 
laser-atom physics, although there is no explicit evidence for their key assumption that the 
ionized electron can be treated by the classical mechanics directly 6|. 

Bohmian mechanics (BM) 7, 8, 9| also called quantum trajectory method is another 
version of quantum theory, in which trajectory concept is used to describe the motion of 
particles with the Bohm-Newton equation. It has been successfully used to study some 
fundamental quantum phenomena such as the two-slit experiment and tunneling 



Recently, Orio 



et al 



12 



13|. 



14| used BM to simulate the electron transport in mesoscopic systems. 



Sanz et al jl5| applied this theory to the atom surface physics. Makowski et al derived 



some central potentials and two-dimensional noncentra 



investigate the exact classical limit of quantum mechanics 
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potentials [17[ from BM to 
. BM has also been regarded 



as a resultful approach to studying chaos [19|, due to its description of trajectory for the 
quantum system. It has also been successfully applied to the intense laser-atom physics to 
study the dynamics of above-threshold ionization |2l| and high-order harmonic generation 



22] by authors. 

The only difference between the Bohm-Newton equation and Newton equation is that 
there is an extra term in the Bohm-Newton equation, called quantum potential. When 
quantum potential is negligible, the Bohm-Newton equation will reduce to the standard 
Newton equation and then the motion of particles can be described by classical mechanics. 

In this paper, BM is introduced to the intense laser-atom physics. We first obtain the 
electron quantum trajectories of an atomic ensemble from the Bohm-Newton equation. Next 
we study how the quantum potentials change as the electrons are driven away from the parent 
ion by the laser field. We find that quantum potentials trend to be smaller and then become 
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completely negligible soon after the electrons were ionized. In this case, the Bohm-Newton 
equation reduces to Newton equation, and then the motion of electrons can be described 
by classical equation of motion. Thus, our results present a direct positive evidence for the 
previous semiclassical trajectory methods in intense laser-atom physics where the motion 
of ionized electrons is treated by classical mechanics directly, while quantum mechanics is 
needed before the ionization. On the other hand, our numerical results clearly show how far 



the electron from the core as it was ionized, while the initial position of the ionized e 

f 

is usually chosen to be (the position of the core) in the semiclassical approaches [2 



ectron 



II. QUANTUM TRAJECTORIES FORMALISM OF BOHMIAN MECHANICS 

BM 0, Isl is derived from a subtle transformation of the time-dependent Schrodinger 
equation. Firstly, a wavefunction can be written in the polar form ip{r,t) = i?(r, t)e*'^*^'''*'*/'^, 
where R and S are real functions. Secondly, the wavefunction is inserted into the time- 
dependent Schrodinger equation. The real part of the resulting equation has the form 

f + (Z^ + v-i!lX!^ = o (1) 

dt 2m 2m R ^ ^ 



and the imaginary part is 



^ + Vipv) = 0, (2) 



where p{r,t) = R^{Y,t)^ v = VS{r,t)/m, and V is the ordinary potential. Equation ([Tj) is 
similar to the classical Hamilton- Jacobi equation, except it has an extra term, This 
term, denoted by Q{r,t) in this paper, is usually called the quantum potential. Equation 
([2]) looks like the classical continuity equation. So a Bohm-Newton equation of motion for 
a Bohmian particle can be constructed: 

md^r/dt^ = -V{V + Q) (3) 

from the standpoint of classical mechanics. The motion of particle is determined by the 
ordinary potential V and the quantum potential Q which plays a crucial role for the ap- 
pearance of quantum phenomena. Obviously, when quantum potential Q ^ in equation 
([3]), the Bohm-Newton equation will reduce to the standard Newton equation and then the 
motion of particle can be described by classical mechanics. 
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FIG. 1: (color online), (a) Graphic description of the direction of the linearly polarized laser field 
(E||z) at frequency to and the coordinate of the Hydrogen system; (b) Laser field profile in the first 
ten periods, where the laser period T = lOvr in our work. 



In fact, the trajectorjM) 
instead of equation ([3] 



if_particles can be evolved from a much simpler equation of motion 

^ = V5(r,t)/m. (4) 

To obtain the trajectory, we need to know the phase S{r, t) and the initial positions of 
particles. According to BM, the initial distribution of particles in an ensemble is given by 
|^(r,0)p. 
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III. INTENSE LASER- ATOM INTERACTION 



The system we consider in this paper is a Hydrogen atom in intense laser field. The 
Schrodinger equation for the system can be written as (atomic units are used throughout) 

t^^^ = mr) + Hjiv,tmr,t). (5) 

Here -f^o(r) is the field-free Hydrogen atom Hamiltonian and Hj{r, t) is the intense laser-atom 
interaction: HQ{r) = — |^ + ^ ~ }:^Hj(v,t) = — r ■ E(t) = —zE{t), where the laser is the 
linearly polarized field (E||z) and E{t) is the laser field profile. Due to the linearly polarized 
laser field, magnetic quantum number m of the atom is a good quantum number, so that 
the problem of solving the time-dependent Schrodinger equation here can be simplified into 
a two-dimensional problem (see figure [l](a)). In our study we use the grid method and the 
second-order split-operator technique to numerically solve the time-dependent Schrodinger 
equation which has been detailedly introduced by Tong and Chu 23|. The laser field 
profile is 

^^^^ ^ I Eosm\§) sin(^t), < t < 3T 
I Eq sin(u;t), t > 3T, 

where T = 2it/uj, Eq and u are the electric field amplitude and angular frequency, respec- 
tively (see figure [U^b)). Here we take Eq = 0.1 au and u = 0.2 au The initial state ilj{r,0) 
of the system is in the ground state of the field-free Hydrogen atom. 



IV. RESULT 



A. Electron trajectory from the Bohm-Newton equation 

After numerically obtaining the time-dependent wavefunction ip{r, t), and then the phase 
S{r,t), we can get the trajectory of an electron by integrating equation (j4]) with its initial 
value of position ro (see appendix A). In this way, we can gain an ensemble of electron 
trajectories according to the initial distribution of electrons |'?/'(r, 0)|^, where 'ip{r,0) is the 
ground state of the field-free Hydrogen atom. Explicitly, we present four trajectories of 
electron with the corresponding initial positions in figureO Let's take one electron trajectory 
as an example with ro = (1.0 au, 7r/30) in the polar coordinate. The path of the electron has 
the following character in the spatial coordinate: The motion is irregular when the electron 
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FIG. 2: (color online). The electron trajectories in the spatial coordinate with the corresponding 
initial positions tq. The initial positions ro in the polar coordinate are from top to bottom: (1.0 
au, 7r/30), (2.0 au, vr/SO), (1.0 au, vr/lO), and (1.0 au, 27r/15), respectively. 

is near the parent ion; after the electron has been driven far away from the parent ion by 
the intense laser field, it travels in a straight line with a little oscillation. The reason why 
the motion of the electron keeps oscillating as the electron is far away from the parent ion is 
due to the existence of the laser field. The other trajectories with different initial positions 
have similar character as shown in figure [2l 

In figure [3l we show the projections of electron trajectory in z and x directions, respec- 
tively, as functions of time. In figure [3](a), two curves z{t) and x{t) are obtained from the 
electron trajectory with ro = (1.0 au, vr/SO). At the beginning, both z{t) and x{t) are near 
the core. After the time t ~ 130 au, z{t) becomes a regular wave curve and x{t) seems 
to be a straight line. Similarly, we have obtained the projections of three other electron 
trajectories with their initial positions tq (see figures 3(b)-(d)). All of them have the similar 
character as that in figure [3](a). 
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FIG. 3: (color online). Projections of electron trajectories in z and x directions as functions of 
time. The corresponding initial positions of electrons are, from (a) to (d), (1.0 au, vr/SO), (2.0 
au, 7r/30), (1.0 au, vr/lO), and (1.0 au, 27r/15), respectively. It is obvious that the motion of the 
electron in x direction is uniform while its motion in z direction has undulations as the electron is 
far from the core. The reason for the z-undulation is due to the existence of laser field (E||z). 

B. Time-dependent quantum potential in intense laser-atom physics 

In Section [TTl we have discussed that quantum potential (5(r, t) together with ordinary 
potential \^(r, t) dominates the motion of particle in BM. And quantum potential plays a 
crucial role for the appearance of quantum phenomenon. If the value of Q{r,t) smoothly 
trends to be negligible (Q — > 0), equation ([3]) will reduce to Newton equation and then the 
motion of electron can be described by classical mechanics. 

Now we study the change of quantum potential for an atomic electron as it is driven away 
by the laser field. After obtaining the time- dependent wavefunction ip{r,t) and the electron 
trajectories r(t), we calculate the changes of quantum potential Q{r(t)) (see appendix B). 
In figure m we show the time- dependent values of quantum potential Q{r{t)) and ordinary 
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FIG. 4: (color online). The time-dependent value of quantum potential and ordinary potential of 
the electrons with the corresponding initial positions, (1.0 au, 7r/30), (2.0 au,7r/30), (1.0 au, vr/lO), 
and (1.0 au, 27r/15), respectively. 

potential V{r(t)) = — ^ — zE{t) with the corresponding initial positions fq. Here we take 
the electron with tq = (1.0 au, vr/SO) as an example (see figure 4(a)). The value of quantum 
potential is comparable with that of ordinary potential in the early period. After some time, 
e.g., after 100 au about, quantum potential smoothly becomes smaller and then it trends to 
be negligible, particularly after the time t ~ 133 au The changes of quantum potential of 
three other electrons are presented in figures 4(b)- (d), and they have the similar character 
as that in figure ll](a). In this way, our numerical results present an explicit picture of the 
reduction from the Bohm-Newton equation to Newton equation. Consequently, the Bohmian 
trajectory method can give both the quantum and classical descriptions of the motion of an 
atomic electron in intense laser-atom physics. 
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FIG. 5: The time-dependent total energies of an electron (including kinetic energy, Coulomb 
potential and quantum potential) with the initial positions, (1.0 au, vr/SO), (2.0 au, vr/SO), (1.0 
au, 7r/10), and (1.0 au, 27r/15), respectively. 

C. Classical description after the electron has been ionized 

As we have shown above, the Bohm- Newton equation of electron is approaching to Newton 
equation in intense laser-atom physics, as the electron was driven far away from the parent 
ion. In the previous semiclassical approaches of intense laser-atom physics, the ionized 



electron motion is usually treated 



used before the ionization 



w the classical dynamics, while quantum mechanics is 



2J]. Here, we want to verify whether the motion of 
an electron can be treated by classical dynamics or not, i.e. whether the quantum potential 
becomes negligible or not, as it has been ionized. 

In figure [5], we show the total energies of an electron (including kinetic energy. Coulomb 
potential and quantum potential) with different initial positions. When the total energy 
becomes positive, the electron will be ionized. Here we take the electron with ro = (1.0 au, 
7r/30) as an example (see figure 5(a)). At the beginning, the value of the total energy is less 
than zero, i.e. the electron is bounded by the core. After the time t ^ 103 au, the total 
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energy becomes positive but it then turns to the negative at the time t ~ 120 au The reason 
why the total energy of the ionized electron turns to the negative again may be that the 
electron with low energy is not far enough from the core and it is recaptured by the parent 
ion with a photon emitted when the laser field reverses. The total energy becomes positive 
again and higher than before at t ~ 133 au and thereafter, i.e. the electron can never be 
recaptured again by the parent ion. On the other hand, in figure |l](a), we found the value of 
quantum potential of the electron trends to be trivial soon after the time t ~ 133 au (The 
corresponding distance between the ionized electron and the core is 16.2 au about; see figure 
3). Similarly, in figures 5(b)-(d) and 4(b)-(d), we can find when the electron is ionized, the 
corresponding quantum potential trends to be negligible. Therefore, our numerical solutions 
show that the quantum potential is negligible as the ionized electron was driven far away 
from the core and thus the motion of ionized electron is dominated by classical mechanics. In 
this way, our work may present a direct positive evidence for the semi-empirical assumption 
in the semiclassical approaches of intense laser-atom physics, such as the three-step model 
[2I and Feynman's path- integral approach in the strong field approximation 24], in which 



the motion of ionized electrons is treated by classical dynamics directly. Furthermore, the 
Bohmian trajectory method can help us to find how far the electron from the core when 
it was ionized, while the initial position of an ionized electron is assumed to be in the 
semiclassical approaches Q|. 

V. DISCUSSION 

As it has been discussed above, quantum potential dominates the quantum behavior 
of a physical system. In the prior work of quasistatic model of strong-field multiphoton 
ionization j^, a dual procedure was considered to describe the multiphoton phenomena in 
intense laser-field physics. First a tunneling model is applied to obtain the ionization rate 
and describe the formation of a sequence of wave packets. The second part of the quasistatic 
procedure uses the classical mechanics to describe the evolution of an electron wave packet. 
It has been pointed out definitely that only the Newton equation cannot give the ionization 



rate correctly, unless a quantum tunneling process was allowed 25| . This is consistent with 
our numerical results in this paper: The quantum potential cannot be ignored before the 
ionization, i.e. quantum mechanics is needed to study the behavior of the electron before 
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the ionization, but the motion of an ionized electron can be described by classical mechanics 
directly . 



VI. REMARK 



It is interesting that Mahmoudi et al [26| have studied free-electron quantum signatures 
in intense laser fields by looking for negativities in the Winger distribution. They concluded 
that the quantum signatures of the free electron get washed away, provided the dipole 
approximation is made. In this way, their work presents another evidence for the treatment 
of ionized electron in the semiclassical approaches of intense laser-atom physics. 
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APPENDIX A: ELECTRON QUANTUM TRAJECTORY 

In this appendix, we show how to obtain the electron quantum trajectory from equation 
(jl]). First, equation (jl]) can be expressed as (atomic units are used) 

where '?/'(r, t) = R{r,t)e^^^^'^^ and Im(/) is the imaginary part of /. In our study '?/'(r, t) is 
gained by numerically solving the time-dependent Schrodinger equation ([5]), using the grid 
method and the second-order split-operator technique, which has been detailedly introduced 



by Tong and Chu [23|. Then we use the Runge-Kutta method to evolve equation ( lAll) to 
obtain the electron quantum trajectory. 

In our numerical procedure, the range of the variable r in the radial direction is confined 
to (0, 150 au). In the grid method, the numbers of grid points are 350 in the r direction and 
51 in the 6 direction. The time step for the evolution is 0.01 au 
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APPENDIX B: TIME-DEPENDENT VALUE OF QUANTUM POTENTIAL 



After obtaining the time-dependent wavefunction ipi'^^t) ^^^id the electron quantum tra- 
jectory r{t) with the corresponding initial position tq, we can numerically obtain the value 
of the quantum potential Q(r(t)) along the trajectory r{t). First, the quantum potential 
Q{r,t) can be written as 

2" 



Q(r,t) 



1 V^i?(r,t) 

2 R{r,t) 



Re 



V^^(r,t) 
^(r,t) 



Im 



V^(r,t) 



(Bl) 



where Re(/) is the real part of /. Secondly, inserting the quantum trajectory r{t) into 
equation (IBip . we can numerically get the value of quantum potential Q{r(t)). 
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